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Abstract 

The conserved Sine-Gordon Equation with nonconserved shot noise is 
used to model homoepitaxial crystal growth. With increasing coverage 
the renormalized pinning potential changes from strong to weak. This is 
interpreted as a transition from layer-by-layer to rough growth. The asso- 
ciated length and time scales are identified, and found to agree with recent 
scaling arguments. A heuristically postulated nonlinear term V 2 (V/i) 2 is 
created under renormalization. 



1 Introduction 

In Molecular Beam Epitaxy (MBE) it is possible to control the amount of de- 
posited matter through oscillations of the surface roughness, which indicate that 
the surface grows layer by layer |], |, |3|. A simple picture represents layer- by- 
layer growth on a high symmetry surface as follows: After deposition out of the 
beam atoms diffuse on the surface until they either meet other diffusing atoms to 
form a stable island, or get incorporated at the edge of a previously nucleated is- 
land. If most atoms deposited on top of an island are assumed to be incorporated 
into its edge by performing a downward hop (which implies that the suppression 
of interlayer transport by Ehrlich-Schwoebel-barriers M is negligible), little nu- 
cleation occurs in the second crystal layer before the first layer is completed. 
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Consequently the surface width at layer completion is nearly zero, after having 
gone through a maximum at half filling of the first layer. As long as the layer- 
by-layer growth mode persists, the surface morphology exhibits oscillations with 
a period given by the monolayer completion time. 

In general this scenario is only transient, and the oscillations are damped. A 
variety of mechanisms contribute to the damping: The surface may be slightly 
miscut [0, |5|, or the average beam intensity may be inhomogeneous ||. How- 
ever, even in the absence of such (experimentally unavoidable) imperfections, the 
stochastic beam fluctuations are sufficient to destroy the temporal coherence of 
spatially separated regions on the surface. Provided the beam noise is the sole 
damping mechanism, it was recently shown that the critical coverage 9 at which 
the oscillations disappear scales with the ratio of the surface diffusion constant 
D s to the deposition flux F as || |7j 

9 ~ (D s /F) s (1) 

with an exponent 

' = "<—* (2) 

where d denotes the surface dimensionality (d = 2 for real surfaces) and the 
exponent 7 characterizes the dependence of the diffusion length (or typical island 

size)£ D onZVF|], 

£ D ~ {D s /F)\ (3) 
Eq.(|l]) may therefore be rewritten as 

9 „ ~ l d (4) 

where 

is the coherence length, an estimate of the size of coherently oscillating regions 



4/(4- d ) (5) 



The theory of Ref. J7| is based on a phenomenological stochastic continuum 
equation for the growing surface. It is not clear a priori that such a continuum 
description would be able to capture the phenomenon of growth oscillations, 
which is distinctly a lattice effect. As a first step towards a more complete 
treatment, in the present work we therefore perturbatively include the lattice 
structure by analyzing the driven, conserved sine-Gordon equation 

d t h = -KA 2 h- \A(Vh) 2 - V Asm — + F + T] (6) 

a± 

for the surface height h(pc, t). In this equation the constant F denotes the average 
deposition flux, while the noise ?7(x, t) models its "shot noise" fluctuations. The 
noise is assumed to be Gaussian with mean zero and correlator 

(r/(x, t)??(x', £')> = 2D 5{t - t')5(x - x'). (7) 
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To motivate the systematic terms on the right hand side of eq.(^|), we note 
that it can be written in the form of a continuity equation 



d t h + V-3 = F + r] (8) 

reflecting the absence of desorption and defect formation under typical MBE 
conditions [0, ITTJl [TT]|. The current is given by Fick's law, J = — Vp, where the 



quantity 

p(x,t) = -KAh- \(Vh) 2 - Vsin— (9) 

a± 

can be interpreted as the spatially varying part of a coarse grained adatom density 
or chemical potential ||11|| . The first term incorporates a generalized Gibbs- 
Thomson effect |12[], according to which the density is enhanced near maxima 
and decreased near minima of the surface, while the second term reflects the 
dependence of the density on the local vicinality [|7], [11], The third term 



models the lattice potential: Adatoms are preferably driven to places where they 
can be incorporated such that the surface remains at integer multiples of the 
vertical lattice constant a±. Only the lowest harmonic of the periodic surface 
potential is kept, since components of periodicity a±/2, a±/3, . . . are less relevant 
(see below and |14}1 ). It is necessary to distinguish the vertical (a±) and horizontal 
(a\\) lattice constants, since they play very different roles in the renormalization 
group calculation. 

The analysis of Ref. was based on eq.([|) with V = 0. Here we show that, 
by explicitly including the lattice potential, the characteristic time and length 
scales (f],^D emerge naturally in the renormalization group (RG) flow equation of 
the potential strength V. Moreover the (Vh) 2 nonlinearity in (Q) is seen to be 
generated under renormalization, thus relieving us from the task of postulating 
its microscopic origin; we may set A = microscopically. A similar scenario is 
valid for the nonconserved Sine-Gordon model [HJ. These results are obtained 



by applying the Nozieres-Gallet RG scheme flTJJ to eq.(J6|). An RG analysis of (y) 
was previously presented by Tang and Nattermann fllEfl , however these authors 
considered separately the cases V = 0, A ^ and V ^ 0, A = and thus were 
not able to address the generation of A from the lattice potential; in addition, 
our analysis includes explicitly the flat initial condition of the surface, which is 
essential for describing transient behavior. 

Since the interpretation of the RG results depends crucially on relating the 
"mesoscopic" coefficient K to the microscopic length scale £d, the next section 
will address this issue within the framework of the linear equation (V = A = 0). It 
turns out that the mere existence of a vertical lattice constant a±_ is sufficient for 
the nontrivial time and length scales to emerge from the continuum theory, 
even if this scale has no dynamical effect (i.e., V = 0) [XT]- The full problem with 



V 7^ is treated in the following sections. After briefly recalling the RG scheme, 
the renormalization of the conserved Sine-Gordon Equation in 2 + 1 dimensions 
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is carried out in Section 3. The RG-flow equations for the parameters in Eq. (||) 
are interpreted in Section 4, and the extension to general dimensionalities and 
general relaxation mechanisms is briefly addressed. Some conclusions are given 
in Section 5. 



2 Length scales in the linear theory 

On the most basic level, the growing surface morphology evolves in response to 
the competition between disordering beam fluctuations and smoothening surface 
diffusion. The simplest continuum theory that incorporates both effects is the 
linearization of (^]), 

d t h = -KA 2 h + F + 7]. (10) 

As was mentioned in Section 1, the coefficient K arises from an expansion of 
the local adatom density in the surface curvature Ah. Under near-equilibrium 
conditions, it would therefore be expected to be given by the product of the 
surface stiffness and the adatom mobility Hl2| , [T~8|1 . However, far from equilibrium 
other processes may contribute to, and in fact dominate K 0, In particular, 



it has been suggested p(J that random island nucleation produces a contribution 

K~F& D , (11) 

but the underlying microscopic mechanism is not known. In the following we show 
how this relation follows from a simple reinterpretation of flKp in the presence of 
a finite vertical lattice constant a±. 

The straightforward solution of ( |T0D [ |TT| ] shows that, starting from a flat 
substrate at time t = 0, after time t surface correlations have developed up to a 
scale 

£{t) « (Kt) 1/4 (12) 

and the surface width grows as 

W(t) « {D/Kf/ 2 i{tf (13) 
in dimensionalities d < 4, where 

C = ^ (14) 



is the roughness exponent of the linear equation ||11|| . 

Together with the average growth rate F the presence of the vertical lattice 
constant induces a fundamental time scale, the monolayer completion time 

tml = a ± /F. (15) 
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Setting t = tml in (0) we obtain a corresponding lateral length scale £(tml), the 
scale on which lateral structure has developed after deposition of one monolayer. 
Clearly it is very natural to identify this scale with the diffusion length £ D , and 
hence the coefficient K in (Tl2f) can be identified as 



K « 4/tml = aJ x F4 (16) 

in accordance with (11). The correlation length flI2|) can then be expressed in 
terms of the coverage 9 = Ft/a± as 

f (t) « ^# 1/4 . (17) 

To see how the scaling laws ([!]) and (||) can be obtained along similar lines, 
note first that for shot noise the noise strength D is proportional to the beam 
intensity F, 

D « cuafF. (18) 

Thus QT3D takes the form 

« a ± (a||/£ D ) d/2 ^ (4 - d)/8 . (19) 

If one now postulates, plausibly 0, that the lattice effects disappear when the 
roughness due to long wavelength fluctuations becomes comparable to a±, the 
characteristic coverage 9 can be defined through W(9) ~ a± and is given precisely 
byeqs. (§ and ©. 

These considerations may be viewed as a zeroth order assessment of lattice 
effects, to be justified by the systematic calculation provided in the remainder 
of the paper. They easily extended to general linear equations with a dynamical 
exponent z ||1 1|| , in which case one finds 



7 (20) 

2 — a 



for z > d. 



3 Renormalization Group Analysis 

3.1 The Renormalization Scheme of Nozieres and Gallet 

Because of the structure of Equation (§) we use an approach which is suitable 
for general forms of the nonlinearity. It was introduced by Nozieres and Gallet 
for the dynamical renormalization of the Sine-Gordon equation to obtain the 
roughening transition |14} . A detailed presentation can be found in their work, 
which we recall briefly. 
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Consider a Langevin equation 



d t h = Ch + M{h) + 7] (21) 

with a linear part Ch, a nonlinear term M{h) and Gaussian noise 77 with mean 
zero and correlator (r/(k, t)r?(k', t')) = 2|k| 2 ^D5(t - t')^(k + k')0(|k| - A). The 
cutoff A = 1/ay is introduced to model the lateral lattice structure which does 
not allow for fluctuations on scales smaller than the horizontal lattice constant 
ay. In the following a\\ will only appear in the cutoff A, hence we can disregard 
the distinction between an and a± and set a± = a. 

Two types of noise can be considered: Either volume conserving noise, which 
corresponds to the case fi = 1 ("diffusion noise" |J), or nonconserving "shot" 
noise fi = 0. Here we focus on the nonconserved contribution, /i = 0, which 
always dominates on scales larger than the diffusion length t D [0, §]. A study of 
the conserved case has been presented in Ref . [^l[ . 



Renormalization of Equation ( pip is performed in the following way: 

• We average over the short wave components 5r] of the noise. In k-space 5r) 
is nonzero only for modes k with (1 — dl)A < |k| < A. Define the averaged 
or coarse grained field h = {h)^. 



• Equation ( pl|) is split in two; one for the coarse grained field 

d t h = Ch + (M(h + Sh)) Sri + fj 

and a second one for the difference Sh = h — h. One now takes an approx- 
imation of (Af(h + 5h))s v in terms of the (hopefully all) relevant operators 
appearing in Equation (|21|). For this one calculates 5h (respectively its 
correlation functions) perturbatively in the nonlinearity M . Since Af is not 
of a simple polynomial form, a Rayleigh-Schrodinger expansion is used - 
the only one feasible one, albeit poorly controlled. 

• Time, lateral and vertical length are rescaled with different exponents: x — > 
(1 — dZ)x, h — > (l — (dl)h and t — > (1 — z dl)t. This causes a corresponding 
rescaling of the terms in ( pip yielding the RG flow equations for C and N ' . 

The detailed application of these steps to Equation (||) is the subject of the next 
section. 

3.2 Application to the Conserved Sine— Gordon Equation 

Consider Equation (||) in a frame moving with the average growth speed F 

8 t h = -KA 2 h - V A sin ^ (h - Ft) + rj. (22) 
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Epitaxial growth starts on an atomically flat surface, so at time t = the initial 
configuration is h = 0. Fluctuations are caused by the noise at later times. This 
will play an important role in the interpretation of our results. 

We neglect the perturbative contribution of the — AA(V/i) 2 -term. As argued 
in the Introduction on a microsopic level it is absent. It is generated by the lattice 
potential and the driving force F to order V 2 . Although it is a relevant operator 
at the linear fixed point, we expect it to be negligible on small and intermediate 
scales, as long as the lattice potential contributes to renormalization of K and 
A. The same is observed for the (V/i) 2 -nonlinearity in the nonconserved Sine- 
Gordon Equation flT5| . 

Averaging over 5r] leads to two coupled equations, where we expand to lowest 
order in the infinitesimal quantity 8h 



dth 



d t 5h 



-K A 2 h - VA 



2vr f- 



sm 



-K A 2 5h - VA 



cos 




(Sh 2 ) 



+ 6r]. 



+ rj 



(23) 



The second equation of (]23|) is solved by a Rayleigh-Schrodinger perturbation 
ansatz 

<5/i (0) (x,t) 



<5/i (1) (x,t) 



J d d x' J dt' G(x - x', t - t') 5rj(x', t') 
2TrV - J d d x' J' dt' G(x - x', t - t') x 
cos — (h(x!,lf)-Flf) ^ (0) (x',t') 



(24) 



a 
A' 



In k-space the linear propagator G is given by exp —K(t — t')|k| 4 . In the sequel 
we will use the notation h', rf, . . . when the argument is the integration variable 
(x',t') and unprimed symbols for quantities at (x, t). The corresponding deriva- 
tives (Laplace operators) are denoted A and A' to mark this difference. For the 
convolution integral J d d x' Jq dt' the shorthand / is used. 

Insertion of (5h 2 ) = (5h^ 2 ) + 2(5h^5h^) + 0(V 2 ) in the first equation of 
(f23|) generates the corrective terms for h to order V 2 . 



(6h 2 ) = (5h^ 2 ) - — [ (A'G) cos — (h' - Ft') C 
a J a \ J 



2tt 



(25) 



where C = (5h^6h^ } is the unperturbed correlator. Its Fourier transform is 



given by e 



-K{t-t')M 



D/(K\k\ 



1-e 

2. 



-2ift|k| 4 



The lowest order term {5h^ ) is a constant which multiplies the sine poten- 
tial, yielding the correction 

4ir 3 D 



SV 



Ka 2 A 2 



-2KAH 



dl V. 



(26) 
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In general dimension 27r/A 2 is replaced by SdA d ~ 4 , where Sd is the surface area 
of the d- dimensional unit sphere. 

To handle the next order contribution due to (Sh^Sh^), 



8tt 3 V 2 



A J S i n ^l(h-Ft\ cos^fh'-Ft'^j (A'G) C 



we split the product sin a cos (3 = (sin(a + /3) + sin(a — (3))/2 and discard the 
term with sin(a + /3). It would create higher harmonics of the lattice potential, 
which are less relevant (see Eq. fl26|)) than the fundamental. In the remaining 
expression we split into terms with h and Ft 



cv' 



A 



sin —Fit' -t) cos — [h-h' ) (A'G) C 
a a 



2tt 



+ cos — F(t' - t) sin 



Expanding terms with h — h' in powers of £ 
symmetries when integrating over <i 2 £ we get 



h'j (A'G) C . 
x' — x and taking care of the 



Att 5 V 2 



2^V 2 



A(Vh) 
A 2 h 



[ sin \—F(t' - t)] (A'G) C £ 2 
J la 



cos 



2tt 



F(t' - 1) 



(A'G) C £ 2 



(27) 



Up to now we have expanded the averaged nonlinearity (J\f(h))s v to order V 2 in 
projections onto the relevant operators. 

The last step consists in rescaling space and time. To examine the behavior 
close to the linear fixed point V = 0, A = we choose its scaling exponents, 
z = 4 and ( = 2 — d/2 = l in two dimensions (see (|HP). Accordingly the growth 
rate F and the lattice constants a± and a» are rescaled as dF/dl = (z — ()F, 
daj_/dl = — (aj_ and dau/dl — — a». 



3.3 Flow equations 

The RG flow will be examined in terms of dimensionless quantities, in which the 
effect of the trivial rescaling has been eliminated. As elementary length scales 
we use the vertical and horizontal lattice constants a± and an, and the basic 
time scale is given by the monolayer completion time tml = Q>±/F. Thus time 
is measured through the coverage 9 = t/r ML , and the coefficients K, V and A 
appearing in (j^) are replaced by the expressions 
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Due to the conserved form of the equation of motion (^j), the noise strength D is 
not renormalized [l6j , and need not be considered further. 

Using the scaling parameter K = exp4Z, the flow equations then take the form 

(i) = -Arc 3 ^ (l - e" 2 ^) U, 

(h) = 9W*,0), (28) 

R 4k ^ = 4^ /(K/M)- 
gsk; /C 

The functions / and g depend on the integrals in Eq. (|27D and are given in the 
appendix. Equations fl28|) are the central result of this paper. The following 
sections are devoted to the discussion of their physical content. 



4 Interpretation of the RG Results 

As in the case of the roughening transition [H|, the renormalization of U de- 
termines the relevance of the lattice on large scales. In the present situation 
the lattice potential always decreases under renormalization (eq.(|28|(i))). Indeed, 
since in the absence of the lattice potential the roughness exponent (|14| ) ( > 
in two dimensions, the lattice becomes irrelevant asymptotically; a roughening 
transition is possible only if ( = 0, such as for eq.flBI) subject to conserved noise 



2j| . Nevertheless on finite length or time scales U remains finite, and its depen- 
dence on 6 and k may be used to describe the transition from (lattice-dominated) 
layer-by- layer growth to rough, continuous growth. It can be shown that for small 
U the amplitude of the surface width oscillations is proportional to U ; thus the 
value of U on a given time or length scale is a direct measure of the observable 
signatures of layer-by-layer growth. 

In the following we assume that the rate of particle deposition is small com- 
pared to the diffusion rate, which is true for typical MBE conditions and implies 
that the dimensionless stiffness parameter /C^>1. The renormalization of K, due 
to the lattice potential, as expressed by eq.(^8|(ii)), can then be disregarded, and 
/C becomes a constant. This decouples the flow equation for the lattice potential 
U and allows for a straightforward solution, which can be used to extract the 
damping time 9 and the coherence length i. In Section EOI the generation of the 



A(V/i) nonlinearity, as described by eq.(p8|(iii)), will be discussed. 
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4.1 Damping time 

The asymptotic (k — > oo) value of the lattice potential U is obtained from 
eq.(|28|(i)) as 

t/oo = e- /(e) f/ (29) 



with 

7r 3 r°° dn 



m = lcl ^ (1 " e ' 2K9/K) w ( 3 °) 

in the relevant regime O ^> 1. Thus writing 7(0) = \Jo /6 the characteristic 
coverage is obtained as 

^=^~(Wa||) 4 , (31) 

07T' 

in agreement with ([|) for d = 2. Moreover, since in this case the surface width 
of the linear theory is proportional to (9/9) l / A (see Section |2|), we see that the 
decay of the lattice potential is of the form 

Uoo/Uq « exp[-CW 2 (#)] (32) 

where C > is a constant. This behavior has been found to describe the decay 
of oscillation amplitudes in layer-by- layer growth in numerical simulations, and 
can be derived assuming a discrete probability distribution of the heights tak- 
ing at each possible height the value of the corresponding continuous Gaussian 
distribution [^2] . 



It is instructive to extend these results to other dimensionalities and linear 
relaxation mechanisms with a general dynamic exponent z. For general z, the 
natural scaling variable is k = expzl. Then the only qualitative change of the 
expressions discussed above is that the algebraic part of the integrand in ( |30"D be- 
comes k 2< '/ z ~ 1 instead of n~ l l 2 . To see this, note that the only quantity changing 
under rescaling in the time-independent part of the correction 5V/V in eq.(PS[) 
is the vertical lattice spacing a_i_ ~ k~^^ z . Thus, for £ < the integral (^) con- 
verges even when 9 — > oo, implying that the surface remains smooth. For £ > 
the integral becomes 

1(9) ~ Kr\K6)^ ,z ~ W\ (33) 

showing that (|32"1) remains valid in the general case. Writing 1(9) = (9/9) 2 ^ z , 
the characteristic coverage is of the order 

9~ K z/K - 1 = )C d/{z - d \ (34) 

where in the last step the scaling relation z = d + 2( for linear growth equations 
[pT| has been used. This becomes equivalent to (p0|) by noting that (|Tj]) is replaced 
by K « £ z d /tml for general z. 
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4.2 Coherence length 

For 9 — > oo the solution of the flow equation (p8|(i)) at finite k is 

U(k) = exp[-(27r 3 //C)(v^ - 1)]U Q . (35) 

If the renormalization is stopped at a lateral scale L» = k 1 / 4 ^, the remaining 
value of the lattice potential may therefore be written, for L\\/a\\ 3> 1, as 

U(L ll )/U = e- {L ^ /i)2 (36) 

with 

i=a n [£ { ~*k, (37) 

in agreement with the expression ([|) for the coherence length. If the system is 
smaller than I the lattice potential remains relevant, the surface remains smooth 
and growth oscillations will be present for all times |7|]. 

As in the previous section, these considerations can be generalized to arbitrary 
z and d. Then fl36J) becomes U(L\\)/U = exp[-(L||/£)^] with 

I ~ /C 1 / 2 ? = K} /{z - d) ~ (Va||) 2/(2 " d) (38) 
in accordance with ( pD| ) and Ref.0. 



4.3 Generation of the conserved KPZ nonlinearity 

We now focus on the flow equation (p8|(iii)), which describes the generation of 



the nonlinear term of the conserved Kardar-Parisi-Zhang p3| (CKPZ) equation, 
A(V/i) 2 , through the interplay of the lattice potential V, the growth rate F and 
the effective stiffness K. We consider the stationary regime 9 — > oo, and again 
assume that /C is large, so that its renormalization can be neglected. Setting 



L = A = at the microscopic scale, the solution of the flow equation (|28|(iii)) 
then reads 

L{k) = ^Lj^xx- 1 / 2 exp[-(47r 3 //C)(v^- l)]f(K/x), (39) 



where the solution (|35|) for the flow of U has been used. For k — > oo, eq.(|39[) 
tends to a finite limit, which for /C — ► oo has the simple form 

Loo = 2vr 2 [/ 2 . (40) 

The numerical solution of the full set ( p8|) of coupled flow equations shows that 
the limiting value (|40|) is attained for /C > 10 8 , corresponding to a diffusion length 
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Id = lOOaii. In the range 10 4 < K, < 10 8 the nonlinear coupling Loo is positive 
and increases with increasing /C, while for smaller values of /C it is negative. 

In terms of the bare parameters of the original equation, the relation ( f40|) 
implies that, on large scales and for large diffusion lengths, the CKPZ coefficient 
A is positive and of the form 

x ~^ ± < 41 > 

It is instructive to compare this to a heuristic estimate of A, based on Burton- 
Cabrera- Frank (BCF) theory |7], 11]. According to BCF [23 , the adatom density 



on a terrace can be computed by solving a steady state diffusion equation with 
sinks at the surface steps. The density decreases with decreasing step distance 
or increasing tilt, and becomes independent of the tilt when the step distance 
is of the order of the diffusion length ip. The coefficient A of the leading order 
expansion @ around a singular surface is then positive^ and given by 

Abcf « (42) 

To identify eqs.(|4l|) and (^) we would need to require that the bare pin- 
ning potential Vq depends on the diffusion length as Vq ~ F£ D , and hence the 
dimensionless potential strength is 

Uq « (W^i) 2 = v^> 1, (43) 

which is clearly inconsistent with our perturbative treatment of the potential. 
Thus, the expressions ( f41~|) and (f42|) are not equivalent, but rather correspond to 
different limiting situations: Our calculation is an expansion for small Uq and 
fixed (large) /C, while the BCF picture assumes perfect crystal planes, which 
would be formally represented by taking Uq — > oo at fixed /C. Writing = 
Loo(U ,IC) we have shown that = Uq 0(/C) for Uq — > 0, where is an 
increasing function of /C, and the BCF argument indicates that ~ JC for 
Uq — > oo. Clearly the latter regime is not accessible by our method. Nevertheless 
it is remarkable that the CKPZ coefficient emerges from the RG calculation with 
the correct sign and the correct qualitative dependence on the diffusion length. 



5 Conclusions 

In this work we have derived renormalization group flow equations for the con- 
served Sine-Gordon Equation with nonconserved shot noise. They were used to 
model the crossover in homoepitaxy from layer-by-layer growth to rough growth. 

1 For a vicinal surface growing in the step flow mode one can show that A < 0, see |L3|. 
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The crossover can be quantitatively characterized by a characteristic layer coher- 
ence length t and an associated coverage 9, which is a measure of the number of 
growth oscillations that can be observed under optimal growth conditions (that 
is, in the absence of miscut or beam inhomogeneity) . The dependence of these 
length and time scales on growth parameters is in agreement with dimensional 
analysis and numerical simulations 0. 

Our approach also predicts the presence of a nonlinear term of the form 
A(V/i) 2 on large scales, which was suggested previously on heuristic grounds 
|| [10| . By power counting it is seen to be relevant in the long time limit and it 
has nontrivial effects on the scaling behavior |25| . 



Two extensions of this work seem to be possible within the conserved Sine- 
Gordon ansatz: First, renormalization of a tilted surface (as performed for the 
original Sine-Gordon model in Ref.[|l4|) should clarify the influence of a small 
miscut on the damping of growth oscillations |2|. Second, and more ambitiously, 
the implementation of an Ehrlich-Schwoebel-effect Q may provide a systematic 
approach to computing the surface current induced by step edge barriers || 
|TT| p6[ , and thus contribute to understanding the transition from layer-by- layer 
growth to a coarsening mound morphology [JIT], p7[ . 

Acknowledgements: We thank H. Kallabis for helpful discussions. This 
work was supported by Deutsche Forschungsgemeinschaft within SFB 237 Un- 
ordnung und grosse Fluktuationen. 



A Appendix: Scaling Functions 

The functions used in the flow equations (|28|,i) and (|28|,ii) are 

i _ p -2K.6/k 

9 {K/k,9) = 8j b 



f(K/K,d) 



[(/C/k) 2 + 7T 2 ] 3 

IC/k(-(1C/k) 4 - Ati 2 (K/k) 2 + 5tt 4 ) 
+ cos(27r#) + A sin(27r#)^ e " 2 ° /K 

i _ „-2H/k 

167T 6 

[(/C/ K ) 2 + 7T 2 ] 3 

7r(^-(/C/fi:) 4 -87r 2 (/C/fi:) 2 + 7r 4 ) 



+ \~ A cos(2vr#) +B sin(27r£ 

with the polynomials 

A = 4tt# 2 (/C/k) 6 + (7t-87t^ 2 )(/C/k) 4 
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+ 87T 3 #(/C//€) 3 + (87T 3 - AlT b d 2 ){]C/Kf 

+8tt 5 9JC/k-tt 5 
B = -46 2 (JC/k) 7 + 46(JC/k) 6 + (1-87i 2 9 2 )(JC/kY j 

+ 16tt^(/C//€) 4 + 4(tt 2 - tt 4 # 2 )(/C/k) 3 + 12^6(JC/k) 2 

— 57T 4 /C//t. 
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Figure 1 

Numerical solution of L(k^ oo) as defined in Equation ( p9[) using the full flow 
equations (|28|). The potential strength is [To = 0.1; for large values of /C > 10 8 
the value 2tt 2 Uq P0| ) is attained. 
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